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The  influence  of  finite  ion  Larmor  radius  and  equilibrium 
self-electric  fields  on  the  ion  resonance  instability 
in  a nonneutral  plasma  column  is  examined,  and  a closed 
algebraic  dispersion  relation  for  the  complex  eigenfrequency  u> 
is  obtained.  It  is  shown  that  finite  ion  Larmor  radius 
effects  can  have  a strong  stabilizing  influence  for  azimuthal 
mode  numbers  £.>2,  particularly  when  the  equilibrium  self- 
electric  field  is  sufficiently  weak. 
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I.  INTRODUCTION 


One  of  the  most  basic  Instabilities  that  characterizes  a nonneutral 

plasma^  with  both  ion  and  electron  components  is  the  ion  resonance 
2-5 

instability.  In  cylindrical  geometry,  the  ion  resonance  instability 

2 

can  be  described  as  a two-rotating-stream  instability  in  which  the 

relative  rotation  between  electrons  and  ions  is  produced  by  the 

equilibrium  self-electric  field  E^(r),gr.  Previous  analyses  of  this 

2-4 

instability  have  been  based  on  macroscopic  cold-fluid  models 
in  which  the  ion  and  electron  motion  is  assumed  to  be  laminar.  In 
this  paper,  we  Investigate  the  influence  of  finite  ion  Larmor  radius 
and  equilibrium  self-electric  field  effects  on  the  ion  resonance 
instability  in  a nonneutral  plasma  column  with  rectangular  electron 
and  ion  density  profiles.  The  analysis  is  carried  out  within  the 


framework  of  a hybrid  Vlasov-fluid  model.  The  electrons  are  described 

as  a macroscopic,  cold  (T^O)  fluid  Immersed  in  a uniform  axial 

magnetic  field  B^e^  On  the  other  hand,  to  allow  for  the  possibility 

of  large  ion  orbits  with  thermal  Larmor  radius  comparable  to  the 

radius  of  the  plasma  column,  we  adopt  a fully  kinetic  model  for  the 

ions  in  which  the  ions  are  described  by  the  Vlasov  equation.  An 

important  conclusion  of  the  present  analysis  is  that  the  growth  rate 

of  the  ion  resonance  instability  exhibits  a very  sensitive  dependence 

on  r../R  , u>2  /w2  and  f.  (Here  rT.  is  the  characteristic  thermal  ion 
Li  p pe  ce  Li 

Larmor  radius,  R is  the  radius  of  the  plasma  column,  u _ is  the 
P pe 

electron  plasma  frequency,  w is  the  electron  cyclotron  frequency, 

C6 

and  f*n?/n^  is  the  fractional  charge  neutralization  provided  by  the 
i e 

positive  ions.) 

For  azlmuthally  symmetric  electron  equilibria  (3/36*0  and  3/3z*0) 
characterized  by  electron  density  n°(r)  and  mean  velocity  ^g(jj)“V®0(r)^e, 
equilibrium  force  balance  on  an  electron  fluid  element  in  the  radial 
direction  can  be  expressed  as  -meV^  2(r)/r—eE®(r)"eVg0(r)BQ/c,  or 


, 


t-xfct--' v JU*.  • '■  * 
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2 0 
equivalently  w*-u)ce where  w^^eBg/m^c,  o>e(r)*Ve0(r)/r,  and  i*>g(r) 

is  the  angular  J^*J£q  frequency  defined  by  Wg^-cE^/rBQ.  For  purposes 

of  analytic  simplification  in  the  stability  analysis,  we  specialize 

to  the  case  of  a sharp-boundary  equilibrium  in  which  the  electrons 

have  a rectangular  density  profile,  i.e.. 


i i 
! 


I! 


\\ 

1 1 

■ 


(n_»const.  , 0<r<R  , 

B®(r)  ■ P P U> 

(o  , Vr<Rc  » 

where  r-R  is  the  radial  location  of  a grounded  conducting  wall, 
c 

For  the  ion  equilibrium,  we  make  the  particular  choice  of  distribution 
function  f^  that  also  gives  a rectangular  density  profile,  i.e., 

fJ-(fnomi/2TT)6(H1-u)iP0-fi)G(vz)  , (2) 

* 2 2 0, 

where  a^,  f,  and  are  constants,  H1»mi(vr+v0)/2+e<))  (r)  is  the 

perpendicular  energy,  P0*m^(rv0+r^uc^/2)  is  the  canonical  angular 

momentum,  and  G(v  ) is  the  parallel  velocity  distribution 
z 

with  normalization  dv  G(v  ) - 1.  After  some  simple 

z z 

i~C°  6 

algebraic  manipulations,  it  is  straightforward  to  show  that 

the  electron  and  ion  density  profiles  precisely  overlap  with  a common 

- 2 2 2 2 
radius  Rp  provided  ^-m^n  Rp/2,  where  hk  is  ion  mass,  0 “^E^ci^i'^i^'ci’ 

and  w -eBrt/m. c.  In  the  case  of  a sharp-boundary  equilibrium,  the 
Cl  U 1 

rotation  frequency  can  be  expressed  as  a)E“2irnQec(l-f ) /Bq.^ 

II.  ELECTROSTATIC  STABILITY  PROPERTIES 
In  the  stability  analysis,  flute  perturbations  with  3/3z*0  are 
considered.  For  perturbations  with  azimuthal  harmonic  number  £, 
a perturbed  quantity  6 ( x , t ) can  be  expressed  as  6i<((^,t)»^  (r)exp(i(i9-u>t) } , 
where  a>  is  the  complex  eigenfrequency.  In  the  electrostatic  approxima- 
tion with  6^(^,t)*-^6«t>(^,t) , the  linearized  Vlasov-fluid  and  Poisson 
equations  can  be  expressed  as^ 


i 
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)V  -(-CD  +2(D  )V  .-(e/m  )(3/3r)$(r)  , 

e er  ce  e eo  e 

-idu-icD  )V  +[-cd  +(l/r)(3/3r)(r2io  ) ]V  -(e/m  )(U/r)*(r)  , 

6 6u  CB  c c 

-i(cD-£<De)ne+(l/r)O/3r)(rt:°Ver)+(lA/r)n°Vee-0  , (3) 

n f°  - 

f,(r,v)-(e/m.Vi)(3/3V1)f:^(r)+i(cD-£cD1)  dx4>  (r' )exp{  iX  (0  '-6)-1wt}  ] 

1^1  * J —re. 

[ (1/r)  (3  /3  r)  (r3  /3  r)-«.  2/r2  ]$  (r)— 4*  e(J d3vf  i-ng)  , 


where  the  perturbation  amplitudes  V , VeQ,  ng  and  refer  to  radial 
velocity,  azimuthal  velocity,  electron  density,  and  ion  distribution 
function,  respectively.  Making  use  of  Eqs.  (1)  and  (3),  and  3n®(r)/3r= 
-n.6(r-R  ),  we  obtain  the  eigenvalue  equation 

0 p 


- 2 

- 0)  —2d  cd 

lPte)  £ _E£  5 ,r_R  ) 

r w-lw  2 (r  V 
e v r 


(4) 


- ^1*  dv*  t ii(r)+("-t“i)ii  • 

2 2222 

where  v (r)»(co  -2u  ) -(cd-JIcd  ) ■const.,  cd  (r)»4irn..e  /m  ■const.,  for 

e ce  e e pe  0 e 

2 

0<r<Rp,  and-(Dpe(r)*0  otherwise,  and  orbit  integral  I is  defined  in 
Eq.  (5). 

It  is  evident  that  the  perturbed  electron  contribution  to  the 
right-hand  side  of  Eq.  (4)  [the  term  proportional  to  6(r-R in  Eq. 
(4) ] is  equal  to  zero  except  at  the  surface  of  the  plasma  column 
(r*Rp) . Moreover,  it  can  be  shown  that  Eq.  (4)  supports  a class  of 
solutions  in  which  the  perturbed  ion  density  [the  term  proportional 
to  jd^v  Vj^3f®/3VA. . . in  Eq.  (4)]  is  also  equal  to  zero  except  at 

r*Rp.  It  then  follows  from  the  linearized  Poisson  equation  (4) 

» „ 1 

that  the  electrostatic  potential  #(r)  has  the  simple  form  $(r)*Ar 


~rT~. .L ' — — - 
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inside  the  plasma  column  (0<r<Rp),  where  A is  a c^astant.  Making 

use  of  the  ion  trajectories  in  equilibrium  fields,  it  is  readily 

. 6 

shown  that  the  required  orbit  integral  can  be  expressed  as 

fO  + 

X-i [$ (r)  / ^ ] J dTexp(-iu>T)[(u)i-u.i)exp(iujiT) 

(5) 

-(u^ui^expdoj'r)  , 


where  u>ci[l±(l+4o>E/a>ci)1/2]/2.  An  important  feature  of  Eq.  (5) 
is  that  the  orbit  integral  I is  independent  of  perpendicular  energy 
®jV2/2.  This  is  a consequence  of  the  particularly  simple  form  of 
$(r)  within  the  plasma  column.  The  right-hand  side  of  Eq.  (4)  is  equal  to 
zero  except  at  the  surface  of  the  plasma  column  (r*Rp) . Moreover,  the  eigen- 
function <fr(r)  satisfies  the  vacuum  Poisson  equation,  r *(3/3r) [r3$/3r]  - 


d2/r2)^(r)»0,  except  at  r-R  . Therefore,  the  solution  to  Eq.  (M 

can  be  expressed  as  $(r)-Ar*  for  0< r<Rp  and  $(r)-Ar  (1-Rc  /r  )/(l-Rc  /Rp  ) 

for  R <r<R  . Note  that  4>  (r)  is  continuous  at  r-R  . 
p - c r 

The  dispersion  relation  that  determines  the  complex  eigenfrequency 
u,  is  obtained  by  multiplying  Eq.  (A)  by  r and  integrating  from  Rp(l-e) 


to  R (1+e)  with  e-KJ. . This  gives 
P + 

1 


1-(R  /*c)21  2(u)-£ue)  [(uj-iwe)-(u)ce-2ije)  ] 2iv‘ 


_E£_ 


u2  R2 


(6) 


where  “p,'4'' “o'2/V  “c.5*V"«c-  and 


r4(co)— 1 + 


‘V’t 


f ! 


w-lui , 


“i^i 


l 

m!(i-m)f  - , + I + 

ra-0  w-iu>^-m(u>^-u>^) 


(7) 


A striking  feature  of  the  present  analysis  is  the  fact  that  the  required  orbit 
integral  i [Eq.  (5)]  can  be  evaluated  in  closed  form  [Eq.  (7)]  for  general 

A a 2 2 

values  of  the  parameters  r_./R  and  (2tu  /id  )(l-f).  Moreover,  the  resulting 
r Li  p pe  ce 

eigenvalue  equation  (4)  for  the  perturbed  electrostatic  potential  $(r)  can  be 

solved  exactly  to  give  a closed  algebraic  dispersion  relation  [Eq.  (6)] 

for  the  complex  eigenf requency  <d.  As  expected,  in  the  limit  where 

r ,/R  +0,  Eq.  (6)  reduces  to  the  familiar  cold-fluid  dispersion  relation 
Li  p n 

2 3 

previously  discussed  in  the  literature.  * 

Equation  (4)  has  been  solved  numerically  for  the  complex  eigen- 

frequency  oj*wr+iY  for  a wide  variety  of  plasma  parameters.^  For 

present  purposes,  we  assume  that  the  electron  fluid  is  rotating  in  the 

1 — 1/2 

slow  equilibrium  mode  with  u)  *id  [1-(1-4<d_/co  ) ]/ 2,  while  the  mean 

n e e ce  F- 


E ce 


equilibrium  motion  of  an  ion  fluid  element  corresponds  to  the  slow  rotation 

-2  2 1/2 

velocity  defined  by  (*)^*-wc^[l-(H-4ajg/u)c^-4r^^/Rp)  ]/2.  The  present 

analysis  is  restricted  to  nonneutral  proton-electron  plasmas  (m^/me»1836> , 

and  the  growth  rate  and  real  frequency  are  measured  in  units  of 

1/2 

the  lower-hybrid  frequency,  w =(ui  id  ) . Moreover,  we  assume 

LH  C6  Cl 

2 2 

that  Rp/Rc»0.5.  Stability  boundaries  in  the  parameter  space  (f,  <*>pe/^ce) 

are  illustrated  in  Figs.  1 and  2.  In  Fig.  1,  the  solid  curves  correspond 

to  the  stability  boundaries  (y*0)  obtained  from  Eq.  (6)  for  Ju'Y0’ 

and  several  values  of  azimuthal  harmonic  number  1.  For  a given 

*2  2 

value  of  £.,  the  region  of  (f,w  /<d  ) parameter  space  above  the  curve 

pe  ce 

corresponds  to  instability  (y>0) , whereas  the  region  of  parameter 

space  below  the  curve  corresponds  to  stability  (ym0) . In  the  forbidden  zone 

in  Fig.  1,  equilibrium  is  not  allowed  since  the  magnetic  restoring 


force  on  an  electron  fluid  element  is  weaker  than  the  repulsive  space- 

*2  2 

charge  force  [2id  (l-f)>io  ]. 
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In  Fig.  2,  the  stability  boundaries  are  illustrated  for  r.  ./R  -1. 

LI  p 

A a 2 2 

Evidently,  for  such  large  value  of  r.  ./R  , the  region  of  (f,u>  /u>  ) 

Li  p pe  ce 

parameter  space  corresponding  to  allowed  equilibria  becomes 
increasingly  limited  by  the  equilibrium  constraint^  that  the  pressure 
gradient  force  on  an  ion  fluid  element  be  weaker  than  the  confining 
electric  and  magnetic  forces  (see  the  uppermost  forbidden  zone  in 
Fig.  2). 

The  dependence  of  stability  properties  on  fractional  charge 

neutralization  f is  illustrated  in  Fig.  3,  where  the  normalized 

*2  2 

growth  rate  is  plotted  versus  f for  “pg /a,ce“0.002  [Figs.  3(a) 

~2  2 

and  3(b)]  and  5 [Fig.  3(c)],  and  several  values  of  mode 

number  l.  We  also  assume  rT  ./R  «0  in  Fig.  3(a),  and  rT  ./R  -0.5  in 

Li  p Li  p 

Figs.  3(b)  and  3(c).  Several  features  are  noteworthy  in  Fig.  3. 

First,  the  number  of  unstable  modes  increases  rapidly  as  f is 
increased.  Second,  for  £>2,  the  instability  growth  rate  is  significantly 
reduced  by  finite  ion  Larmor  radius  effects,  particularly  when  the 
equilibrium  self  electric  field  is  sufficiently  weak.  [For  example, 
compare  Figs.  3(a)  and  3(b)  with  f approaching  unity].  Third,  the 
instability  growth  rate  increases  substantially  with  increasing  plasma 
density  [Figs.  3(a)  and  3(b)]. 

III.  SUMMARY  AND  CONCLUSIONS 

Several  important  conclusions  follow  from  the  present  analysis. 

For  example,  the  growth  rate  of  the  ion  resonance  instability  exhibits 

a very  sensitive  dependence  on  rT./R  , i2  /u>2  and  f.  Moreover, 

Li  p pe  ce 

finite  ion  Larmor  radius  effects  can  have  a strong  stabilizing 

influence  for  mode  numbers  £>2  [see,  for  example.  Figs.  3(a)  and 

3(b)],  particularly  when  the  equilibrium  self-electric  field  is  weak 
«2  2 

(<Dpe/a>ce<<l  or  f close  to  unity).  For  the  fundamental  mode  (£*1), 
however,  stability  properties  are  identical  to  those  calculated 


L I 


from  a macroscopic  two-fluid  model,  and  the  growth  rate  is  unaffected 

by  the  value  of  r,  ./R  . The  detailed  dependence  of  normalized 
Li  p 

*2  2 6 
real  frequency  on  r^/R^,  “pe^ce  an<*  ^ 8^-ven  elsewhere. 
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FIGURE  CAPTIONS 

Fig.  1 Stability  boundaries  [Eq.  (6)]  in  the  parameter  space 

(f,i2  /oi2  ) for  m,/m  -1836,  R /R  -0.5,  r ./R  -0,  and 
pe  ce  l e pc  l-i  p 

several  values  of  £. 

Fig.  2 Stability  boundaries  [Eq.  (6)]  in  the  parameter  space 

(f ,a>2  /w2  ) for  m,/m  -1836,  R /R  -0.5,  r../R  -1.0,  and 
pe  ce  i e p c Ll  p 

several  values  of  Si. 

Fig.  3 Plots  of  versus  f [Eq.  (6)]  for  1^/11^-1836,  Rp/Rc-0.5. 


